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^ Abstract. Integrable multi-component lattice equations of the Boussinesq family have 
been known for some time. Recently some new equations of this type were found using the 
Consistency- Around-the-Cube approach. Here we investigate one of these models, B-2, and 
in particular the consequences of a nonzero deformation parameter 6o > 0^ which allows 
special kinds of solitons in the parameter range — 6o/3 < k < bo. 
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1 Introduction 

> 

^ The lattice Boussinesq equation (IBSQ) reads [1] 

^ y — xx-\-z = 0, y — xx + z = 0, (1.1a) 

y — xx-\-z = 0, y — xx + z — {)^ (1.1b) 



y — XX -\- z — — — - — U, (1-lcj 

where we have used the standard shorthand notation, e.g., x — Xn+i^m^ ^ — ^n,m+i5 and where p 
and q are parameters associated with the n and m directions, respectively. Equations (1.1a), 
(1.1b) are defined on the edges of the elementary square of the Cartesian lattice, while (1.1c) is 
defined on the square itself. These equations can be naturally extended to a third dimension 
and the extension is Consistent- Around-the-Cube (CAC) [1]. (For further information about 
CAC see, e.g., [2, 3, 4].) 

The deformation called B-2 in [5] is given by 

y — bo{x — x) — XX z — — — - = 0. (1-1^0 

The set (1.1a), (1.1b), (l.lc^) also has the CAC property [5]. By changing the sign of x, 

if necessary, we may assume without loss of generality that bo > 0. This will be used in the 

following. 



"^This paper is a contribution to the Proceedings of the Conference "Symmetries and IntegrabiUty of 
Difference Equations (SIDE-9)" (June 14-18, 2010, Varna, Bulgaria). The full collection is available at 
http://www.emis.de/journals/SIGMA/SIDE-9.html 
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In order to understand the physical meaning of the new parameter 60 in we construct its 

sohton solutions. The method is essentially the same as the one given in [6], where we considered 
the standard IBSQ. We repeat here the salient parts of that construction, concentrating on the 
differences. 



2 Construction of the one-soliton solution 

In the construction we use the consistency cube of Fig. 1 in an essential way. Here F — (x, z)^ 
and r stand for the direction parameters of the three directions. In practice consistency 




Figure 1. The CAC-cube. 

means that the equations on the sides provide a Backlund transformation between the bottom 
and the top equations. In the present case the side equations are 

- - - ^ z—^ 

y ^ XX X ^ ~ X ^ i^._Lci) 

tXj 00 00 00 

z = boix — x) + XX — y + — — (2.1b) 
z = 6o(x — x) + XX — y + — — (2.1c) 

along with some of their shifts. 

2.1 The background solution 

First one has to construct the background solution. This is operationally the same as in [6], 
but the parameter 60 will induce some changes. We follow the fixed-point idea [7] which means 
that the variables on the bottom and top square are the same. Thus, omitting the bar-shift 
from (2.1), the equations to solve are 

y — X z ^ z z — X ^x X ) , z z — x (^x x ) , 

-^7/-^ X ~ p - r , _ q - r 
z — bo[x — X) + XX — y + — , z — bo[x — X) + XX — y + — . 

00 00 00 00 

The solution is given by 

xo = an + bm + ci, (2.2a) 
= 1x1 + I (a^n + b^m + C2) + C3, (2.2b) 
yo = \xl-\ (a^n + b^m + C2) - C3, (2.2c) 

where we have introduced new parameters a, 6, which are related to q by 

- boa^ =p-r, b^ - b^b^ = q - r, (2.3) 
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and ci, C2, C3 are arbitrary constants (in [6] a different sign convention was used). Here the new 
parameter 60 makes the correspondence shghtly more involved. 

Since the parameters a, b seem to be the natural ones for this equation we may write (l.lc^) 
using them obtaining 



a 



3 



y — XX -\- z — — — — 60 

X — X 



(x — x) — 



62 



X — X 



0. 



Now we observe that the bo term is in fact related to lattice potential KdV equation (IpKdV) 
given by 

(x — x){x — x) — — b^ 

In other words, (l.lc^) can be interpreted as a combination of IBSQ (1.1c) and IpKdV. 

2.2 The one-soliton solution 

Next, in order to construct the one-soliton solution (ISS) we again use the CAC-cube, now 
with variables x, ^, z at the bottom square corresponding to the background solution and the 
bar-shifted variables x, z on the top square to the ISS. Thus, using parametrization (2.3), 
equations to solve are now 

-— - Z — Z ^ Z — ^ 

y — XX — z^ X — — — X — — — (2.4aj 

X — X X — X 

z = 6o(^ - x) + XX - y + ^-^^^ (2.4b) 

X — X 

z 6o(x - x) + XX - y H ^— ^= (2.4c) 

where k is the soliton parameter. 

Following the method given in [6] we expand the ISS as 

(x, y, z) = (xo + X, + Y, zo + Z), (2.5) 

where (xQ^yQ^zo) is the bar-shifted {xo,yo, zq), i.e., 

Xo = an + bm + A: + ci, (2.6a) 
zo = ^xl + l (a^n + b^m + A:^ + C2) + C3, (2.6b) 
Vo = 1^0 - 1 («^^ + + A:^ + C2) - C3. (2.6c) 

With these definitions we find from (2.4a) that Y = xqX. Thus we only need to solve 
for X, Z, for which we have from (2.4) 

- _ -xqX + Z - _ -{zo + yo + boXo)X + (60 + xo)Z 



X 



X-{a-ky X-{a-k) 
-XqX + Z - -(zo + yo + boXo)X + {bo + xo)Z 



x-{b-ky x-{b-k) 

This system can be linearized by taking (X, Z) = (^, j), because then we can write it as 
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where 



xo —1 

N = lzo + yo + boxo -bo - xq 

-1 a-kj 



xo -1 

M = lzo + yo + boxo -bo - xo 
-1 6- 



These matrices satisfy the integrabihty condition NM = MN. An important observation is 
that the M matrices can both be diagonahzed using the matrix Q defined by 



Q{n,m) = 



( xo{n,m) +bo - UJ2 
xo{n,m) + bo- uJi 

-{xo{n,m) + bo- k) 







k{Sk - 2bo) ) 
M = Q{n, m + l)-^D{b)Q{n, m), 



\ k{3k - 2bo) 
because then we have 

N = Q{n + 1, m)-^D{a)Q{n, m) 

where 

/a-L02 

D{a) = a-LOi 

\ a-k^ 

The entries u!i{j^ k) appearing in the above equations are the roots of 

- - bo{W^ - k'^) =0. 



(2.7) 



This always has the root loq = k, and if 6o = then ooi = kQ, i = 1,2, where C 7^ 1 is the 
cubic root of unity. The key observation is that changing bo changes the values of the other 
roots u)i{k), and although the coi, i = 1,2 can still be complex conjugates of each other although 
not of magnitude k, we also have the novel possibility that all Ui can be real. In fact 



ujj{k) 



1 



{bo -k) + {-iy^{bo-k)ibo + 3k) 



-n-. 2 

So coi are real when —bo/S<k<bo- Using the above formulae we can construct 

^{n,m) = Q(n,m)-iD(a)"D(6rQ(0,0)*(0,0), 
from which we find 



(2.8) 



f^j2(a-u,{k)nb-co,{k)rpl, 

1^=0 
2 

g = -k){a- u:.{k)T{b - Mk)rpl, 

2 

Here, instead of qqq, hoQ, /oo, we have introduced new constants p'^, which are defined by 

_ f _ goo(ci + bo- k) - hpo _ ffoo(ci + fep - o^i^) - hpo _ 

Po-Joo ' (_l).(^^_^2)(a;.-fc)' ' 



(2.9a) 
(2.9b) 
(2.9c) 



Soliton Taxonomy for a Modification of the Lattice Boussinesq Equation 



5 



Using (2.5), (2.6) we can recover the ISS as 

PO + Pi + P2 

z =zo + ^o \ \ \ \ \ , (^.lUbj 

PO + Pi + P2 PO + Pi + P2 

y^^'^yo + x,^-^P^±^^P^±^. (2.10c) 

PO + Pi + P2 

Here xq, yo? were defined in (2.2) and 

p,(n, m) = (a - a;,(A:))^ (6 - c^,(A:))^ p^, = 0, 1, 2. (2.11) 

The pi, correspond to the plane-wave factors (PWF) in the continuous case and we use the same 
name here. 

3 The different types of solitons 
3.1 Generic properties 

The sohton (2.10), when expressed in variables x, ^, z, has a linearly or quadratically growing 
background part (2.2). In order to show the soliton behavior we ignore this part and only discuss 
the part 

U,^k^l^^^P^^^^P^^^^P\ (3.1) 

/ PO + Pi + P2 

From (3.1) we can immediately see that the soliton is a kink or anti-kink and that in the 
asymptotic region it will level off to some a;^, depending on the term that dominates in the 
particular asymptotic region of n, m. In comparing the terms i and j it is useful to determine 
the line where they are equal, it is given by 



nlog 



a — ujj 



+ mlog 



h - Ui 



b-UJj 



0. 



Then the i term will dominate in the growing n region with respect to this line, if \a — uji\ > 
\a — etc. Furthermore one has to verify that the i = j term is dominating over the k term, 
and therefore only a half-line is relevant. In the generic case one needs to draw three such 
half-lines, which results in three different asymptotic regions. Examples of these will be given 
below. 

3.2 The behavior of sohtons for bo = 

Let us first consider the classical case for which bo = 0. Then cj^ = kC where ( = - ^(1 + i^/3). 
For simplicity let us assume that < /c < a, 6. We have 

\a — uji^l = \/ a'^ -\- ak k'^ for z/ = 1, 2, 
from which it follows that \pi \ = \p2\. Furthermore let us define 

a + k/2 



^a{k) = arccos 



\/o? + ak -\r k'^ 



so that a — ujy ^ \/ + ak + /c^e^^^^^^. Since uol = UJ2 we can have real solitons if we choose 
(P2)* — Pi — without loss of generality we can also take = 1- 
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Figure 2. Typical oscillation for the 60 = case. On the left a 3D view of U of (3.1), on the right 
a section at m = 0. The parameters have values 60 = 0, a = 6 = 2, /c = 0.5, = 1. 



We can now write the soliton part of (3.1) in the 60 = case as 
A^S^ + 2acos(C- f^) 



U = k- 



A^B"^ + 2a cos ^ 



(3.2) 



where 



A 



a — k 



B 



b-k 



^ = n(pa{k) + m(pb{k) + /3. 



Since < A: < a, 6 we have < A, S < 1. Thus A^B'^ +oc as n, m ^ — oc and therefore 
in that direction U ^ k. On the other hand, A^B^ ^ as n, m ^ +00 and the behavior 
in that direction is oscillatory since the cosine terms dominate. In the continuous case there 
would be some values of the independent variables for which the denominator actually vanishes. 
In the discrete case we can only say that for large enough n, m we can choose their particular 
values so that ^ gets arbitrarily close to {N + l/2)7r. And then due to the extra — 27r/3 in the 
numerator, U can take arbitrarily large (positive or negative) values. Thus in the 60 = case 
all single-soliton solutions are effectively singular asymptotically. This is illustrated in Fig. 2. 

3.3 The behavior of sohtons when k < — 6o/3 or 60 < (&o > 0) 

The solitons behave essentially the same way as for 60 = also in the 60 7^ case, provided 
that k > bo OT k < — 6o/3, i.e., {k — 6o)(6o + 3A:) > 0. The roots cj^, z = 1,2 are still complex 
conjugates, with |c<j^p = k{k — bo), and \pi\ = \p2\- The detailed expressions for variables A, B, 
(fa(k), used in (3.2), will just have a different form. Let us assume 60 < < a, 6, then we have 



UJ2 



bo-k- iy^{k-bo){bo + 3k) 



60 - A: + iy^{k-bo){bo + 3k) 

a — (jJi\ = \/a2 + a/c + A:^ — 6o(<^ + k) 
a + k/2-bo/2 



for i = 1,2, 



^a(A:) 



arccos 



V^a2 + afc + A:2 - bo{a + k) ^ 
For the plane- wave factors we have rations 
a — k _ b — k 



A 



B 



\a — (jJi\ 



\b - cjil ' 



A' 



1 + 



a(6o -k) + k{bo - 2a) 
o? ^ak^k'^ -bo{a + k) 



<1, 
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or more generally 



uJo{k) 



(a - fc)3 



(> 0). 



One example of such a soliton is given in Fig. 3. 
3.4 The behavior of sohtons when —bo/3 < k < bo 

If —bo/ 3 < k < bo the situation changes completely, because the roots ujjy{k) are now all real, 
and different in magnitude 



k-y^{bo-k){bo-\-3k) 



0J2 



bo-k+ y^{bo-k){bo + 3k) 



Now different terms of U can dominate depending on the value of k. 
For the elements appearing in the PWF's We find 



- Ibo < k <0 : 
< k < Ibo \ 
Ibo < k <bo : 



k < Wi < W2, 
wi < k < W2, 
Wi < W2 < k, 



a — k > a — wi > a — W2, 
a — wi > a — k > a — W2, 
a — wi > a — W2 > a — k. 



(3.3a) 
(3.3b) 
(3.3c) 



The relationships of the middle column above can be seen if we plot k, wi, W2 together, see 
Fig. 4. 

There are now a variety of soliton types depending on the values of soliton parameters k, uji 
with respect to the equation parameters a, b. All the PWF's are positive, if a, 6 > max(A:, it;2), 
and if a,6 < miii{k^wi) the negative sign can be canceled from all terms. In these cases there 
will be no oscillations. If however a and/or b falls between min(/c,it;i) and max(/c,it;2) then 
there will be oscillations and furthermore there can be a greater variety in dominant asymptotic 
behavior. We will next discuss some of the possible cases. 



3.4.1 a, 5 > max(fc, k;2) or a, 6 < m\xi{k^w\) 

If a, 6 > max(A:,it;2) all PWF are positive. The term that dominates for large positive values 
of n, m is (a — k)'^ (b — k)'^ for (3.3a) and {a — wi)'^ {b — wi)'^ in all other cases. For large negative 
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-0.2 0.0 2 0.4 0.6 0.8 1.0 



Figure 4. The roots of (2.7) as functions of /c for 60 = 1: W2 in green, wi in red, and ujq = k in blue. 








Figure 5. When a^b > max(/c, (^2) the sohton looks like a simple anti-kink, bo = 1^ k = —0.2, a = 2, 
6 = 3, pO = 1. 

values for both n, m the dominating term is (a — k)'^(b — k)'^ for (3.3c) and (a — W2)^{h — W2)^ 
in all other cases. For other n, m directions we must compare values more closely. Furthermore, 
choosing a relatively large will make the corresponding term dominate in some finite region, 
although not asymptotically. 

In the example of Fig. 5 the a, b values are so large in comparison with k that the solution 
looks like a simple anti-kink. Similarly for large negative values of a, b we get a simple kink. 

If a, b are only moderately larger than 60 the double soliton nature of the kink starts to 
appear. This is illustrated in Fig. 6. A kink is shown in Fig. 7. 

3.4.2 u?! < a, 5 < ujq^u:2 

For this choice the PWF's po and p2 have alternating signs as functions of n, m, but pi > 0. 
This means that at the line where \po\ = pi or \p2\ = pi there will be unbounded oscillations 
in the asymptotic regions, while the behavior is smooth around the line |po| — \p2\- One such 
case is given in Fig. 8, for which the parameters are 60 = 1, A: = 0.5, a = 0.2, b = 0.3, p^ = 1. 
It then follows that cuq — 0.5, ui — —0.309, 002 — 0.809 and a — cjq = —0.3, b — uq = —0.2, 
a — cji = 0.509, 6 — cji = 0.609, a — UJ2 = —0.609, b — 002 = —0.509. The oscillatory transition 
lines are at n = —0.475m for |po| — I Pi I and n = m for |p2| = |pi|, while the transition line 
\po\ — \P2\ at n = —0.758m is smooth. Another example is illustrated in Fig. 9, the opening 
angle to the lowest level is now smaller. 
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u: 




Figure 6. A 3D view and sections for m = —70, —50, —30, —10, 10 of an antikink with 6o = 1, /c = |, 
a = 1, 6 = 2, = 1. All a — uji^ h — uoi are positive. 



u 




Figure 7. A 3D view and sections for m = -70, -50, -30, -10, 10 of a kink with 6o = 1, /c = -0.25, 
a = —0.3, h = —1.5, = 1. All a — uji, h — uOi are negative. 

3.4.3 < a < ujQ < b < u:2 

The behavior is still more diverse when a < cjq < This is illustrated in Fig. 10, where the 
parameters are h{) — 1^ k — 0.5, a — —0.18, h — 0.63. The parameters bo^k^Ui are the same as 
in Section 3.4.2 but since now b > uq the signs are different: a — cjq = —0.68, b — cjq = 0.131, 
a- (jji = 0.129, 6 - cji = 0.939, a - U2 = -0.989, b - 002 = -0.179. Now all the transition 
lines are oscillatory, the line |po| = \pi\ oscillates as a function of n, the lines |po| — \p2\ and 
\pi\ — \pi\ as functions of m. 



3.4.4 ujQ = u;2 

Finally we have an example for cjq — uj2. This leads to a pure kink or antikink solution, because 
Po — P2- This happens, e.g., for 60 = l,fc = 2/3. Since a — ojq < but a — cji > we have 
oscillations on the transition region. This is illustrated in Fig. 11. 
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Figure 9. 3D view with cross section at m = —25. The parameters are bo = 1^ k = a = 0.2, b = 0.12, 
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Figure 10. Soliton solution with parameter values bo = 1, /c = 0.5, a = —0.18, 6 = 0. 63, = 1. Above: 
3D view with a section at m = 10, below: two adjacent sections at m = —14, —13. The levels are, from 
left, 1, and 2. The 0-2 transition line is smooth in n for even m and oscillatory in n for odd m. 



4 Bilinearization, Casoratians and AT-soliton solutions 

Above we have discussed only the ISS. Although the figures look complicated one can obtain 
the multi-soliton solutions using a procedure similar to [6] . We only list the main results 
By the dependent variable transformation 

9 9 h 9 s . 

X^Xq-j, z^Zq-Xoj + j, y = yQ-Xoj + j, (4.1) 

we can bilinearize the B-2 lattice consisting of iiBl=(l.la), iiB2=(l.lb), HB3={l.lc^) as 
HBl = HB2 = 



HB3 = 



ff ff 

BsB4 + (a - b)ffB4 + [a2 + a6 + 6^ _ bo{a + b)]ffBs 



ix-x)ffff 

where we have also used the parametrization (2.3), and where the bilinear equations are 

Bi = f{h + ag) - g{g + af) + /? = 0, (4.2a) 
B2 = f{h + bg) - g{g + bf) + fs^O, (4.2b) 
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Figure 11. A kink with oscihations. Here bo = 1^ k = 2/3^ a = 0.2, 6 = 0.12 and therefore po = p2 = 
(— l)^+^|p2|, while pi > 0. Top: 3D view and a section at m = —3, bottom: a section with n + m = 1 
on the left and n + m = on the right. 

Bs = f9-f9+ {a - h){ff- ff) = 0, (4.2c) 

= [a^ + a5 + 52 _ ho{a + h)] (// - //) 

+ (a + 6 - ho)(fg - ff) + }s + fh-gl^O. (4.2d) 

The set of bilinear equations (4.2) admits A^-soliton solutions in the following Casoratian form, 

/=|]V^|, g^\N^2,Nl h^\N^2,N + 5 = |]V^, TV - 1, 7V|, (4.3) 

composed of '0 = ('^1, '02, • • • , '^A^)"^ with 
2 

^,(n,m,/) = ^^52(-a;3(fc,)y(a-a;,(fc,)r(6-a;,(fc,)r, (4.4) 

where cJo(fc) = ooi{k) for z = 1, 2 are defined in (2.8) and g^^^ are constants. Here a A^th order 
Casoratian of the column vector 

^lj{n, m, /) = ('^Ai(n, m, /), ^lj2{n, m, /),... , ^iv(n, m, l))^ 

w.r.t. shift variable /, is defined by 

C'n,m(^; {^i}) = |'^/^(n,m,/i),^(n,m,/2), • • . , '^(n, m, /at)]. 
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Usually we write this in the shorthand notation [8] \h,l2, • • • Jn], and in particular for con- 
secutive sequences we use M = 0, 1, . . . , M. 

The proof for (4. 2d) is similar to the one in [6]. In fact, the Casoratian column vector 
defined in (4.4) satisfies the shift relation 



'0 — '0 = a'0, — — h^l)^ 

T'^ = '0 + 6o'0, 



(4.5a) 
(4.5b) 



where V is some N x N matrix. We note that the relation (4.5a) is the same as the one in [6], 
and it is sufficient for proving the first three equations in (4.2), so here we do not consider them 
further. The relation (4.5b) is different from the one in [6]. Using it one can generate an explicit 
form for the identity (Tr(r)/)/ = (Tr(r)/)/ (c.f.. Appendix of [6]), 



TV - 5, TV - 3, TV - 2, TV - 1, ^(TV - 2)| - |TV - 4, TV - 2, TV, ^(TV - 2)| 
+ |TV^, TV + 1, ^(TV - 2) I + + 9-bf + bob^'^s - h) 

- b^-^f [|TV^, TV - 3, TV - 2, TV - 1, ^(TV - 2) | - |TV^, TV - 2, TV, ^(TV - 2) 
+ |TV^, TV + 1, ^(TV - 2) I + a^+V + g - af + b^a^-^{s - h)] = 0, 

by which one can then verify (4. 2d). 



5 Conclusions 

In this paper we have discussed the various soliton solutions of the recently discovered [5] defor- 
mation (l.lc^) of the lattice Boussinesq equation (1.1). It is well known that for Boussinesq-type 
equations the soliton solution is formed by a combination of three plane waves, but in the stan- 
dard case with 6o = the parts merge and produce a kink with oscillations. For 6o 7^ the 
three components are more independent and the solutions look like resonating two-kink solu- 
tions with three different asymptotic levels. Depending on the choice of parameters there can 
be oscillations on the transition region between levels. 
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